Doping-dependent study of the periodic Anderson model in three dimensions 
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We study a simple model for /-electron systems, the three-dimensional periodic Anderson model, 
in which localized / states hybridize with neighboring d states. The / states have a strong on-site 
repulsion which suppresses the double occupancy and can lead to the formation of a Mott-Hubbard 
insulator. When the hybridization between the / and d states increases, the effects of these strong 
electron correlations gradually diminish, giving rise to interesting phenomena on the way. We use 
the exact quantum Monte Carlo, approximate diagrammatic fluctuation-exchange approximation, 
and mean-field Hartee-Fock methods to calculate the local moment, entropy, antiferromagnetic 
structure factor, singlet correlator, and internal energy as a function of the f-d hybridization for 
various dopings. Finally, we discuss the relevance of this work to the volume-collapse phenomenon 
experimentally observed in /-electron systems. 

PACS numbers: 71.10.Fd, 71.27.+a, 71.15.Nc, 64.70.Kb 



I. INTRODUCTION 

The periodic Anderson model (PAM) describes the 
qualitative physics of solids in which a set of localized (/) 
orbitals having a strong on-site repulsion hybridizes with 
another set of noninteracting conduction (d) orbitals. 1 
The basic magnetic features encapsulated in the PAM in- 
clude the formation of moments on the localized orbitals 
due to the suppression of double occupancy by the repul- 
sion. The magnetic nature of the ground state is deter- 
mined by the strength of the hybridization between the 
localized and conduction electrons. 2 If the hybridization 
is small, then the magnetic ordering of the localized mo- 
ments takes place through the establishment of a modu- 
lation of the conduction-electron spin density, an indirect 
Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling. 3 For 
strong hybridization the screening of magnetic the mo- 
ments of the correlated electrons by conduction electrons 
(Kondo screening) is expected. While some of these 
general ideas are well accepted, others are still under 
study. 4 For example, Noziere's notion of 'exhaustion', 5 
that only the conduction electrons near the Fermi sur- 
face can participate in singlet formation and hence there 
are not enough of them to form singlets with all the 
local spins, has received considerable recent scrutiny. 6 
Related to this is the effectiveness of the compensating 
cloud around each localized orbital. The size of this 
cloud has been a subject of controversy. Arguments 
based on scaling and renormalization group ideas for 



the single-impurity Kondo point to a cloud with length 
£k ~ vp/Tx, where vp is the Fermi velocity and Tk is 
the Kondo temperature 7 ' 8 . Since Kondo temperatures 
are normally of the order of tens of degrees, this leads to 
a large Kondo cloud of the order of thousands of lattice 
spacings a. On the other hand, similar arguments for 
the single-impurity multichannel Kondo model have led 
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This cloud has never been observed 



experimentally and Knight shift experiments from Boyce 
and Slichter 10 were interpreted to indicate the absence of 
a large cloud. 

The charge excitations of the PAM also show inter- 
esting features. As in the one-band Hubbard model, a 
strong repulsion on the localized orbitals makes them ex- 
hibit the characteristics of a Mott insulator, including a 
suppression of the / density of states and a vanishing /- 
electron compressibility. 11 The density of states can also 
exhibit a Kondo resonance as it evolves to a single broad 
peak at weak coupling. 11 

While many of these general ideas are well accepted, 
much of the quantitative physics of the PAM remains 
open. Only relatively recently has the phase diagram 
of the PAM been studied by non-mean-field approaches 
like the determinant quantum Monte Carlo 12,13 (QMC) 
and dynamical mean- field theory 14 (DMFT). The former 
provides an exact treatment of correlations on finite lat- 
tices but is computationally demanding and therefore has 
so far been restricted in its applications to the particle- 
hole symmetric point, where both the / and d orbitals 
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are half filled. The DMFT approach works explicitly in 
the thermodynamic limit and hence can observe finite- 
temperature phase transitions without recourse to finitc- 
size-scaling techniques, but it involves the approximation 
of ignoring the momentum dependence of the self-energy 

There are several reasons to believe that the removal 
of the restriction of half filling might reveal different 
physics. First, half filling is the most optimal electron 
density for antiferromagnetism. Since a number of the 
potential systems for which the PAM has been suggested 
as an appropriate Hamiltonian do not exhibit such mag- 
netic correlations, it is desirable to study the properties of 
the doped system, for which antiferromagnetism is con- 
siderably suppressed. 

Second, the noninteracting bands have special fea- 
tures which are unique to half filling. For example, in 
the absence of interactions, the canonical choice of a 
momentum-independent hybridization Vk = V, which 
corresponds to a hybridization of the localized orbitals 
with the conduction orbitals on the same site, gives rise 
to a band insulator at half filling. Meanwhile, the choice 
V k = —2V(cosk x + cosky + cosfc z ), which, instead, hy- 
bridizes the localized orbitals with the conduction or- 
bitals on the nearest-neighbor sites, 15 leads to a metallic 
state at half filling; yet, this hybridization Vk vanishes at 
the half-filled Fermi surface and, therefore, is also some- 
what artificial. While the addition of strong correlations 
is expected to modify the physics greatly, it is still pos- 
sible that the properties of the PAM reflect these special 
features of the noninteracting limit at half filling. 

Finally, the particle- hole symmetric point is a very spe- 
cial one for analytic theories. For example second-order 
perturbation theory (SOPT) in the on-site interaction 
becomes exact at the strong- coupling limit, in addition 
to being exact at the weak-coupling limit. This odd re- 
sult is only true at half filling and, indeed, is only true 
for models with a single impurity band hybridizing with 
a single conduction band. Studying the PAM off half 
filling will therefore give valuable insight into the nature 
of analytic approaches to the case of multiband systems, 
where this anomalous success of SOPT does not occur at 
any filling. 

A specific instance where all three issues come into play 
is the Anderson impurity model and the PAM as they are 
applied to the Cerium (Ce) volume collapse. 16-19 While 
recent QMC studies of the half-filled PAM have lent 
greater quantitative understanding of this problem, 20 the 
relevance of this work for the Ce volume-collapse phe- 
nomenon would be much better justified if the same 
physics could be demonstrated in the doped system. 

In this paper we will present a detailed study of 
the magnetic, charge, and thermodynamic properties 
of the PAM for various band fillings. We will make 
use of the exact QMC, approximate diagrammatic 
fluctuation-exchange approximation (FLEX), and mean- 
field Hartree-Fock (HF) methods. We will concentrate on 
the effects of doping of additional electrons into the sys- 
tem for parameter regimes where the number of localized 



(i.e., /) electrons still remains approximately one. This 
is a natural consequence of the Mott-Hubbard gap in the 
PAM and similarly is also the experimental situation in at 
least one application of interest — the Ce volume-collapse 
transition. 

Our key conclusions are the following: (1) The more 
accurate treatment of correlations provided by QMC and 
FLEX shift the HF transition to stronger coupling and 
decrease its sharpness. HF for the PAM gives more ac- 
curate energies at low temperature near strong coupling 
than weak. (2) Antifcrromagnctic correlations remain 
robust for a range of dopings of the conduction band. 
By comparison the one-band Hubbard model loses its 
long-range antifcrromagnctic order more rapidly upon 
doping. 21-23 (3) Anomalies in the singlet correlator, en- 
tropy, and internal energy persist even when the f-d hy- 
bridization has substantial values on the Fermi surface. 
They also persist to doping levels where antiferromag- 
netic correlations have been eliminated. However, their 
sharpness does appear correlated with the presence of 
magnetism. (4) FLEX captures many of the same quali- 
tative features exhibited by QMC in the dependence on 
the interband hybridization of the local moment, anti- 
ferromagnetic correlations, entropy, and internal energy, 
especially at weak coupling, where the HF internal energy 
is not very accurate. At the same time, as is well known, 
FLEX has difficulty reproducing the Mott-Hubbard gap 
in the strong-coupling limit. 24 

II. MODEL AND METHODS 
A. Periodic Anderson model 

We will study the grand-canonical Hamiltonian 

E {M<^U. - tfMjj* + /iUjJ + H.c.} 

+ £ f^2 n f^ - V^2( n dia + Tlfia) 
i, (T i, (J 

+ C/ /E n / i T n /U- (!) 

i 

d\ a and fl create d (conduction) and / (localized) elec- 
trons with spin a on the lattice site i. We consider 
a three-dimensional periodic simple-cubic lattice, where 
(i, j) is a sum over nearest-neighbor bonds. Therefore, 
this term represents the nearest-neighbor hybridization 
integrals (i.e., hoppings) of the d-d (tdd) an d d-f (tfd) 
orbitals. The remaining one-particle terms are £/, which 
measures the /-orbital energy level with respect to that 
of the d orbital (which was chosen to be zero by conven- 
tion), and the chemical potential /i. The charge density 
operator on site i in orbital d with spin a is n<n a = d\ d. 
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and similarly for rif la . Finally, Uf describes an effective 
on-site repulsion between the localized / electrons. 

When the hopping term is Fourier transformed, the 
resulting expression is £ k CT [ek4<Aa + V ^iAaUa + 
ft a d k(T )], with e k = -2t dd (cosk x + cosk y + cosfc z ) and 
Vk = —2tfd(cosk x + cosky + cosfc z ). The particular 
choice of Sf = —Uf/2 and n = results in (rid) = (n/) = 
1 ((rid) = J2<j( n di<?) etc.). This is the so-called "sym- 
metric" PAM, in which both orbitals are half filled re- 
gardless of the choice of all the other parameters and 
temperature. 

We will express all energy units in this paper in terms 
of tdd = 1- We will also set Uf/tdd — 6 through- 
out the paper and vary tfdjtdd for several fillings and 
£f/tdd- For a- and 7-Ce respectively, calculated values 
for the effective- hopping and Coulomb-repulsion param- 
eters in eV are t d d « 0.90 and 0.74, t fd « 0.19 and 0.14, 
Uf w 5.7 and 5.9. 18 Our choices of the model param- 
eters are therefore in rough agreement with the values 
for Cerium. Finer tuning would not be useful given the 
approximations inherent in using a two-band model. 25 

B. Determinant quantum Monte Carlo 

The determinant- QMC method that we employ is 
standard and has been described many times in the 
literature. 26 Its key features are an exact treatment of the 
interactions via a path-integral representation of the par- 
tition function and the decoupling of the interaction via 
the introduction of a Hubbard-Stratonovich field, whose 
possible configurations are then summed over stochas- 
tically. The approach is able to study several hundred 
interacting electrons, much greater than other exact ap- 
proaches like diagonalization, but still representing small 
lattices, especially in three dimensions. This is an espe- 
cially serious consideration in determining the presence 
or absence of long-range order. Local quantities like the 
internal energy and near-neighbor spin correlations have 
much smaller finite-size effects. The simulations reported 
here are on lattices of 64 (4 x 4 x 4) spatial sites with fill- 
ings ranging from (n) = 2.0 to (n) = 3.6 electrons per 
site, where (n) = (rid) + ( n f)- 

Besides finite spatial-lattice size, a crucial bottleneck 
in the determinant QMC is the "sign problem," 27 which, 
roughly speaking for this application, restricts the tem- 
peratures accessible to simulation to those greater than 
T ~ IF/50, where W is the noninteracting bandwidth. 

C. Fluctuation-exchange approximation 

We will also study the PAM using the fluctuation- 
exchange (FLEX) approximation. 28 The FLEX approxi- 
mation can be motivated in various different ways. It can 
be thought of as an extension of the Hartree-Fock theory, 
in which the bare interaction is self-consistently screened 



by the exchange of electron-hole and electron-electron 
(two-particle) fluctuations. Or it can be viewed as a con- 
serving approximation within the framework of Baym- 
Kadanoff's generating-functional formalism. 29 Another 
view is as an extension of the random-phase approxi- 
mation (RPA), in which the two-particle propagators of 
RPA are self-consistently included in the self-energy Yet 
another, but a very elegant one, is that FLEX is the first 
step in the parquet-approximation's 30 two-particle-self- 
consistency scheme. 

FLEX suffers from lack of the two-particle self- 
consistency found in the parquet approximation, which 
rcnormalizes the bare vertices (which are used to cal- 
culate the fluctuation propagators) with vertex correc- 
tions. Although it is possible to formulate FLEX to 
work in spontaneously-broken-symmetry states (with an 
anomalous self-energy), such as a magnet or supercon- 
ductor, FLEX as commonly practiced rarely does so. The 
only such FLEX calculations to this date include those 
for a superconductor. 31 For previous results using FLEX 
for Hubbard-like models, along with comparisons with 
QMC, see Ref. 32. 

The normal-state FLEX solutions of Hubbard-like 
models are also known to give the incorrect strong- 
coupling limit. It is not known if the broken-magnetic- 
symmetry solutions may give the right limit for main 
thermodynamic quantities at strong coupling. 

The FLEX approximation used in this paper incorpo- 
rates all four channels (density, magnetic, singlet, and 
triplet fluctuations). In order to satisfy certain symme- 
tries, all these four channels are needed. In fact we also 
tested FLEX without using the particle-particle fluctu- 
ations, and it turned out that the resulting accuracy 
is lower. We also use the numerical renormalization- 
group method on the Matsubara frequencies to acceler- 
ate evaluation of the FLEX self-energies and fluctuation 
propagators. 33 

III. RESULTS AND DISCUSSION 

In this section we will show how the properties of the 
PAM evolve from strong to weak coupling at several dif- 
ferent fillings. We begin by describing the orbital oc- 
cupations, showing how the charge is transferred with 
the temperature. This helps establish relevant temper- 
ature scales in these simulations. We then show short- 
range magnetic correlations — measures of the local mo- 
ment and singlet formation — before turning to the issue 
of long-range antifcrromagnctism. The section concludes 
with results for the entropy and internal energy and a 
general discussion of the relation of their behavior with 
magnetism. 

As mentioned in Section II A, the PAM Hamiltonian 
in Eq. (1) with £/ = —Uf/2 and fi = has particle-hole 
symmetry — (rid) = (rif) — 1 — he., both the conduction 
and localized orbitals are half filled, regardless of all other 
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parameters. We will compare the properties of the PAM 
at two electron densities with the behavior at half fill- 
ing. The density (n) — 2.2 is only lightly doped away 
from half filling and allows us to examine the rapidity of 
the destruction of the antiferromagnetism. The density 
(n) = 2.6, as we shall see, is large enough that all traces 
of long-range order have been suppressed. In both these 
cases, we have adjusted the / site energy to maintain 
(n f ) ~ 1. 

Finally, we also present a few results for an extreme- 
filling case, in which (rid) ~ 1.8 and (n/) ~ 1.8, where 
the effect of electron-electron correlations are shown to be 
minimal, as is expected since this is equivalent to a very 
low electron density by a particle-hole transformation. 



A. Orbital occupations 

We begin by examining the evolution of the local or- 
bital density with the temperature T in Fig. 1. At high T 
the orbitals are equally filled since energies like Ef, which 
distinguish the orbitals from each other, are washed out 
by large thermal fluctuations. As T is lowered, the oc- 
cupations shift. By T ~ 0.125 ((3 = 1/T ~ 8), much, 
but not all, of the charge transfer from f to d orbitals is 
complete for this value of tfd = 0.8. For smaller tfd the 
occupations reach their asymptotic ground state values 
at higher temperatures owing to the presence of a larger 
Mott gap. Figure 1 also shows the data for the Hartree- 
Fock solutions to the problem, demonstrating that these 
techniques are in qualitative agreement with the exact 
QMC calculations for the temperature dependence of the 
orbital densities, despite the fact that Hartree-Fock pre- 
dicts magnetic order below T ~ Uf/4 (— 1.5 for Uf =6), 
well above any possible QMC transition temperature. 
The FLEX results are in even better qualitative agree- 
ment with the exact QMC results in the entire tempera- 
ture range shown. 



B. Local moment 



Figure 2(a) shows the dependence of the /-orbital mo- 
ment, (m 2 f f) = ((n/ iT - n fil ) 2 ) = (n f ) - 2(n /iT n /u ), on 
the hybridization tfd at low temperature. Here data for 
several dopings are shown. For all three cases, (rriff) 
is fully formed at weak hybridization, and, as expected, 
it is reduced as tfd increases and the system becomes 
less strongly correlated. Even for relatively large tfd, 
however, (m 2 j) is significantly enhanced over its uncor- 
rected value of 

/™2 



jj) = 1/2. For example, the value 
ff) — 0.8 at tfd = 1 and (n) = 2.2 corresponds to a 
double occupation of (^/i-f^/ij.) = 0.10, a value reduced 
by a factor of two-and-a-half from the result in the ab- 
sence of correlations, (n/iT^/u) = ( n fil)( n fil) = 1/4. 
Although the moment is reduced as tfd increases, even 
at large tfd, the moment is well enough formed to allow 
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FIG. 1. The dependence of the / and d orbital occupations 
on the temperature T at tfd = 0.8 for (n) = 2.2. Other 
parameters are tdd = 1 and Uf = 6 (throughout this paper), 
and Sf = —2.785. Also, the lattice sizes in this paper are all 
4x4x4 unless otherwise specified. Results for QMC, FLEX, 
paramagnetic Hartree-Fock (PMHF), and antiferromagnetic 
Hartree-Fock (AFHF) are shown. 



for interesting magnetic ordering to occur. For U/t = 4, 
the half-filled one-band Hubbard model has (m 2 ) = 0.77 
at low T, a value comparable to the smaller of those 
in Fig. 2(a); yet, the one-band Hubbard model still ex- 
hibits clear long-range antiferromagnetic order for this 
U jt. Later in this paper, we will examine antiferromag- 
netism in the PAM. 

In Fig. 2(b) the FLEX results for the local and in- 
stantaneous component of the magnetic susceptibility, 
% m (AR = 0, At = 0), are shown for different lattice 
sizes. 34 This quantity is the equivalent of the local mo- 
ment. As depicted in the figure, the FLEX and QMC 
results are very similar, although FLEX somewhat un- 
derestimates the local moment. The finite-size effects 
are very small, as would be expected for a local quantity. 



C. Singlet correlator 

We will now turn to the issue of the Kondo screen- 
ing of the local moments. Here we probe the formation 
of local Kondo singlets, in which one conduction elec- 
tron and one localized electron in nearest-neighbor sites 
form a bound singlet state. In order to do that we use 
the singlet correlator, which is the nearest-neighbor in- 
terband spin-spin correlation function and is defined by 

c fd = | EH^/i • Sdi), where S fi = ( /,\ ) ff ( & ) 
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FIG. 2. (a) The dependence of the QMC local moment 
on the hybridization tfd at low temperature T = 1/8 
(/3 = 1/T = 8) for (n) = 2.0, 2.2, and 2.6. e/ = -3.0, 
— 2.785, and —2.6, respectively. Other parameters are the 
same as in Fig. 1. (b) The FLEX results for the same param- 
eters for three-different lattice sizes, showing finite-size effects 
are small. The quantity shown is the AR = 0, At — com- 
ponent of the /-orbital magnetic susceptibility, which is the 
same as the QMC local moment shown in (a). 34 The QMC 
and FLEX results are very similar. 



and likewise for Sdj ■ This definition of Cfd is a natural ex- 
tension of the dot product of the local electron spin and 
conduction electron spin on the same lattice site which is 
widely used to measure singlet formation within a PAM 
with on-site fd hybridization. As discussed above, our 
model has an intcrsite fd hybridization which leads to 
singlet formation of the local electron spin with conduc- 
tion electrons on neighboring sites. 20 

In Fig. 3 we show the dependence of Cfd on the hy- 
bridization tfd for different fillings, calculated by QMC. 
There is a qualitative difference between the (n) = 2.6 
case, and the (n) = 2.0 and 2.2 cases. In the latter 
there is a sharp onset in the singlet correlator Cfd near 
tfd = 0.6. This was suggested to be intimately related to 
the Kondo volume-collapse scenario in our previous work 
for (n) = 2.0. 20 This onset moves to smaller tfd and be- 
comes much less sharp at (n) — 2.6. The (n) — 2.0 and 
2.2 cases are almost identical, as will be seen throughout 
this paper. Thus the fact that our choice of hybridization 
vanishes at the (n) = 2.0 Fermi surface as noted earlier — 
but not at the (n) = 2.2 Fermi surface — shows that the 
former artificiality has no impact on the conclusions in 
Rcf. 20. 

Let us call Tk the temperature where the local singlets 
are formed. We expect Ti((tfd) to increase with tfd. 
The behavior of the singlet correlator in Fig. 3 can then 
be explained: the temperature used in the simulations 
(T = 1/8) crosses Tk- As the crossing point depends 
on band filling we observe different onsets, as discussed 
above. The crossing of T and T K also explains why Cf d — > 
for small tfd- 




FIG. 3. The hybridization dependence of the singlet cor- 
relator, which measures the degree of the local antiferromag- 
netism between an / moment and the conduction electron 
spins on the neighboring sites. Parameters are the same as in 
Fig. 2(a). 
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D. Antiferromagnetism 

Another possible phenomenon for the PAM is an- 
tifcrromagnctic ordering of the moments. Fig. 4(a) 
shows QMC results for the equal time antiferromagnetic 
structure factor Sff, the sum of the staggered correla- 
tions of the local / spins across the lattice, Sff(T,q = 

(7r,7r,7r)) = ' Sfi- This structure fac- 

tor is normalized in such a way that if there is only short 
range magnetic order, Sff is of order unity, whereas if 
there are long range correlations Sf / is of order the lattice 
size (and hence diverges in the thermodynamic limit). 

In Fig. 4(b) we show FLEX results for the instanta- 
neous antiferromagnetic susceptibility Xm- Just as the 
instantaneous, local susceptibility is the FLEX equivalent 
of the local moment, Xm is the FLEX analog of the QMC 
structure factor. 34 We will discuss the issue of long-range 
order below. However before doing so, we emphasize the 
key observation that such order for (n) = 2.6 is clearly 
completely suppressed. 

In interpreting these results further, and in particular 
the possible occurrence of long-range order at (n) = 2.0 
and (n) = 2.2, it is important to emphasize that the Neel 
temperature depends on tfd- As in the three-dimensional 
single-band Hubbard model, Tjv is low at small tfd, rises 
to a maximum, and then falls again at large tfd. The 
important point is that a set of simulations at fixed tem- 
perature, like those of Fig. 4, can potentially cut across 
the T/v vs. tfd curve, yielding the impression of the ex- 
istence of AF order only in some intermediate range of 
tfd- With this in mind, we believe the fall-off of AF order 
at small tfd for (n) = 2.0 and 2.2 reflects only that Tm 
has gone below the simulation temperature. Meanwhile, 
the fall-off at large tfd may also reflect this reduction of 
Tn or the competing formation of Kondo singlets. We 
note that in the PAM with on-site hybridization, one 
expects singlet formation to usurp long-range magnetic 
order completely at large tfd- 12 Our FLEX results sug- 
gest the intersite hybridization is different in this regard. 
Only a very detailed study including finite-size scaling 
can sort out these different possibilities. 

To reemphasize, despite all these complications, Fig. 4 
still carries an unambiguous message: the filling (n) = 
2.6 has a very different behavior in its magnetic proper- 
ties: there is no long-range order. We have already seen 
that the singlet correlator still exhibits a growth with t fd 
at (n) = 2.6. Below we will see that, in the same way, the 
(n) = 2.6 entropy and energy still behave rather similarly 
to cases where the densities are close to half filling. 



E. Entropy 

The low-temperature entropy, S(T = Xi ow ), offers a 
means of probing Kondo-singlet formation as well as an- 
tiferromagnetic order. Specifically, the entropy per atom 
of free spin- 1/2 moments is In 2, which can be reduced 
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FIG. 4. (a) The dependence of the antiferromagnetic struc- 
ture factor on the hybridization tfd, as computed with QMC. 
Parameters are the same as in Fig. 2(a). (b) The FLEX re- 
sults for the same parameters for three-different lattice sizes. 
The quantity shown is the Q = Qaf = (7T, it, n), At = com- 
ponent of the /-orbital magnetic susceptibility, which is the 
same quantity as in (a). 34 The QMC and FLEX results are 
very similar. Note the finite-size scaling of the FLEX results. 
That antiferromagnetism persists at large tfd, when singlet 
formation is expected in the case of on-site hopping between 
/ and d orbitals, reflects our near-neighbor choice of f-d hy- 
bridization. Antiferromagnetism is clearly absent for all tfd 
at (n) = 2.6. 
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by the formation of local singlets in cooperation with the 
surrounding valence electrons or by the onset of mag- 
netic order. Since the / and d subsystems decouple at 
small tfd, the free-moment regime is realized by the / 
electrons in this strong-coupling limit for Ti ow above the 
Neel temperature Tn, yet well below some characteristic 
temperature of moment formation. T\ ow = 0.125 is a rea- 
sonable choice since Fig. 2(a) shows fully formed / mo- 
ments at small tfd for this temperature, while Tn — > as 
tfd — * 0. The entropy of the d subsystem may be exactly 
calculated for tf d — 0, and at T low = 0.125 has already 
reached its T = limit: S d (t fd = 0,T = 0) = 0.43322, 
0.29462, and 0.31277 for (n) = 2.0, 2.2, and 2.6, respec- 
tively. These nonzero S d values at T = are finite-size 
effects of the 4x4x4 lattice arising from the degeneracies 
of the d levels at the Fermi surface, which are a nonneg- 
ligible fraction of states in the Brillouin zone for small 
lattices. The / electrons are moved away from the Fermi 
surface by the Mott-Hubbard gap and arc not affected 
by such finite-size effects. We calculate the entropy for 
the fully interacting system from 



,oo l 

S(t fd ,T) = S 00 - J t dT'- 



1 dE{t fd ,T) 



dT' 



(2) 



and 



= 41n4 - (n) ln(n> - (4 - (n» ln(4 - (n», (3) 

using fits to the QMC energies as described 
elsewhere. 20 ^ 35 

In Fig. 5 we show the dependence of the low-T en- 
tropy on the hybridization tf d at T — 0.125. In order 
to make the analysis of the entropy more meaningful, we 
plot S(tf d ,T) — S d (tf d = 0, T = 0) as a function of tf d . 
After this subtraction the entropy curves in Fig. 5 start 
with the value In 2 at tf d — 0, reflecting the free / mo- 
ments. The FLEX results are slightly off, representing 
the weak-coupling nature of the approximation, but it 
is very interesting to see such a weak-coupling diagram- 
matic approximation doing so well at the strong-coupling 
limit, capturing all of the local-moment entropy at the 
lower fillings and 90% at the highest filling. For (n) = 2.0 
the QMC entropy starts to drop around tf d = 0.6; this 
value moves down to a smaller value of tfd = 0.2 for 
(n) = 2.6. One note should be added at this point; as the 
entropy drop for (n) = 2.6 took place at a smaller tf d , we 
didn't pursue calculations for higher t f d values. From Eq. 
2 it is clear that a thin grid of temperatures is needed to 
get the entropy for a given t f d . We were not able to reach 
the needed low temperatures due to "sign problems" , as 
mentioned in section II B. The drop in the entropy is also 
sharper in the (n) = 2.0 case. The (n) = 2.2 results (only 
shown for FLEX) are almost identical to the (n) = 2.0, 
as seen in other quantities as already noted. The FLEX 
and QMC results still don't fully agree as one approaches 
tfd = 2 (weaker coupling) because FLEX becomes truly 
valid only for tf d >3, as we will show in the next section. 
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FIG. 5. The dependence of the entropy on the hybridiza- 
tion tfd- Both the QMC and FLEX results are shown. Pa- 
rameters are the same as in Fig. 2(a), and the tfd = entropy 
of the d system is removed as discussed in the text. The QMC 
and FLEX results agree rather well, especially in the atomic 
{tfd = 0) limit where they produce the In 2 entropy due to 
the formation of local moments. 



Examination of the singlet correlator (Fig. 3) and anti- 
ferromagnetic structure factor (Fig. 4) suggests that the 
sharp drop in the entropy for the (n) — 2.0 and 2.2 cases 
could be related to both the formation of the Kondo sin- 
glets and the onset of the antiferromagnetism. Magnetic 
order, however, is ruled out for the (n) = 2.6 case, where 
there is nevertheless still a significant drop in the entropy 
with increasing tfd- 



F. Internal energy 

In this section we present results for the internal en- 
ergy as a function of t f d at the relatively low temperature 
of T = 0.125 (in one case 0.167). Because the one-body 
part of the Hamiltonian dominates the tf d dependence 
of the internal energy, we report our results here as dif- 
ferences from the Hartree-Fock (HF) — either the para- 
magnetic (PMHF) solution or the broken-symmetry an- 
tiferromagnetic (AFHF) solution — in order to better ex- 
pose the more interesting effects of correlations. This also 
makes sense given that the HF solutions must be upper 
bounds on the correct internal energy at T = 0. More- 
over, we might expect the PMHF solution to be correct 
in the weak-coupling limit [tfd — > oo) and have observed 
previously that the AFHF energy is correct in the strong- 
coupling limit (tfd — > 0) at low temperature. 20 ' 36 With 
these two limits in mind, we first present results as a 
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function of the ratio tfd/{tfd + 1), where values and 1 
correspond to the strong- and weak-coupling limits, re- 
spectively. 

Figure 6 presents our results for the correlation energy 
E — £pmhf as a function of tfd/ (tfd + 1) at four-different 
fillings: (n) = 2.0 (a), 2.2 (b), 2.6 (c), and 3.6 (d). We 
calculate E with three different methods: QMC, dia- 
grammatic FLEX, and mean-field AFHF. In the strong- 
coupling (atomic) limit, AFHF gives the correct energy 
at every filling, agreeing with QMC. 36 Toward the weak- 
coupling limit, FLEX agrees with QMC very well, giving 
not only the correct weak-coupling limit itself but also the 
correct leading-order dependence on 1/t/d- The PMHF 
solution is, in general, correct at the weak-coupling limit 
but as is evident from Fig. 6(c), approaches this limit 
with an incorrect dependence on 

There is an apparent constant offset between the 
PMHF and other curves, including AFHF, at the weak- 
coupling limit in Figs. 6(a) and (b). Here one is at half 
filling or very close to it, and the AFHF solution on small 
lattices is lower in energy at the weak-coupling limit than 
the PMHF. In the thermodynamic limit, these offsets 
vanish, which is also suggested in the FLEX size depen- 
dence tests to be cited shortly. 

In Fig. 6(d) there is no stable AFHF (or other mag- 
netic) solution at any coupling strength. In addition 
E — -Epmhf is very small compared to (a), (b), and (c). 
Indeed correlations are not important in this case of ex- 
treme filling. Note that while the total fillings (n) vary 
in Figs. 6(a), (b), and (c), in each case (n/) ~ 1 and 
the correlation energies are of comparable size. In sharp 
contrast (nj) ~ 1.8 in Fig. 6(d), confirming the general 
understanding that correlations are largest near half fill- 
ing of an interacting orbital and diminish for fillings of 
this orbital away from half filling. 

The FLEX approximation is always paramagnetic if 
one does not include the anomalous-self-energy diagrams. 
Therefore, unlike the mean-field treatments of correla- 
tions, FLEX has no anomalies arising from the solution 
switching from a broken-symmetry (e.g., antiferromag- 
netic) to the uniform (paramagnetic) state as occurs in 
Fig. 6(c). Such HF transitions are absent in Figs. 6(a) 
and (b) since the AFHF solution is always stable. 
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FIG. 6. The difference of the internal energy for various solutions from the PMHF solution, E — Bpmhf, as a function of 
tfd/(tfd + 1), which is and 1 at the strong- and weak-coupling limits respectively. The results for QMC, FLEX, and AFHF 
are shown, (a) (n) = 2.0, (b) (n) — 2.2, (c) (n) — 2.6, (d) (n) — 3.6. In (d) £/ = —5.4, and in (c) (3 = 6. Other parameters 
are the same as in Fig. 2(a). In all cases QMC agrees with AFHF toward the strong-coupling limit and with FLEX toward 
the weak-coupling limit. In the rather extreme filling case of (d), the energy difference E — -Epmhf is two orders of magnitude 
smaller than the other cases, implying that the correlations are negligible and the PMHF solution is quite accurate. In general 
at the weak-coupling limit, PMHF gives the right solution except when the AFHF solution has a lower energy than the PMHF 
as in (a) and (b), which can happen for small finite systems; then, AFHF gives the right solution at the weak- as well as the 
strong-coupling limit. 
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The results in Fig. 6 were for 4x4x4 periodic cells. 
To illustrate the effects of size dependence, Fig. 7 shows 
■Eflex — -Epmhf calculations done for 4 x 4 x 4, 6 x 6 x 6, 
and 8x8x8 periodic cells. As expected the finite-size 
effects on the internal energy are generally small. This 
is because the internal energy is basically a local quan- 
tity. Note, however, the upward shift in Sflex — -Epmhf 
with increased cell size for larger tfd in Figs. 7(a) and 
(b), which correspond to Figs. 6(a) and (b), respectively. 
These shifts are in the direction to remove the constant 
offsets at the weak-coupling limit between the PMHF 
and more rigorous methods in Figs. 6(a) and (b). It is 
also the case that the lowest T = HF solution at the 
weak-coupling limit is PMHF for all fillings in the ther- 
modynamic limit. 
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FIG. 7. The difference of the FLEX internal energy from 
the PMHF internal energy, plotted as a function of the hy- 
bridization tfd for different lattice sizes for the purpose of 
evaluating the finite-size effects, (a) (n) = 2.0, (b) (n) = 2.2, 
(c) (n) — 2.6. Parameters are the same as in Fig. 2. These 
finite-size effects are generally small. 



G. Crossover behavior in the PAM 

We now discuss some of the implications of these re- 
sults for the potential presence of sharp switches in the 
PAM energy which could drive phase transitions in ap- 
propriate systems. 

In the recent work 20 for the PAM at half filling, 
(n) = 2.0, we noted coincident anomalous behavior in 
three distinct low-temperature quantities as tfd was in- 
creased beyond about 0.6: the singlet correlator increased 
abruptly; a Kondo resonance appeared and grew at the 
Fermi level; and the internal energy broke away from 



the AFHF result it had closely followed at smaller tfd- 
We suggested that the abruptness of this change in in- 
ternal energy, evidently associated with the screening of 
the local moments by the conduction electrons, was not 
inconsistent with the physics of the volume-collapse tran- 
sitions in the rare-earth metals. Such a possibility also 
was previously discussed in the context of the Anderson 
impurity model. 16 We now revisit this issue from a dif- 
ferent perspective than given in our earlier work. 20 

Figure 8(a) shows our QMC and PMHF internal en- 
ergies relative to the AFHF results, E — -Eafhf, for all 
three fillings at the low temperature of T = 0.125. Since 
the physical dependence of the hybridization on volume 

is tfd ~ V"~ 2 , we have chosen v = tj^ 2 oc V to label our 
horizontal axis. Note first the AFHF to PMHF transition 
at v = 0.96 (tfd — 1.09) which occurs as v is decreased for 
the filling (n) = 2.6. While there are no HF transitions at 
the two other fillings for the 4x4x4 cell, such transitions 
do occur in the thermodynamic limit at v = 0.87 and 0.74 
(tfd = 1-33 and 1.83) for (n) = 2.6 and 2.2, respectively, 
but there is none exactly at (n) = 2.0. These mean- 
field calculations describe a transition from a "localized" 
AFHF phase, with the / bands Hubbard split away from 
the Fermi level, to an "itinerant" PMHF phase, with the 
merged / bands overlapping the Fermi level, in funda- 
mentally the same way as realistic HF and modified local- 
density methods describe the volume-collapse transitions 
in the rare-earth metals. 18 While having some degree of 
validity, it has been argued in Ref. 18 that the realis- 
tic HF and LDA+U transitions occur too close to the 
itinerant limit for the rare earths, corresponding to too- 
small values of v in the present PAM analog. The QMC 
results in Fig. 8(a) appear to support this suggestion. 37 
While both the AFHF and PMHF values must be upper 
bounds on the exact QMC at T = 0, it is interesting that 
at large v the QMC values are quite close to the AFHF 
upper bound, whereas at smaller v they appear to bend 
and stay well away from the downward-moving PMHF 
bounds. The net effect is a downturn in the QMC ener- 
gies with decreasing v well before the PMHF curves cross 
below the AFHF, as evident for (n) = 2.6 and likely to be 
so in the thermodynamic limit for (n) = 2.2. This is at 
least suggestive of the correlated version of the collapse 
transition, although, a compelling case is hindered by the 
internal-energy contributions omitted by the PAM. 

It should also be emphasized that as v is decreased, 
the onset of the downturns in the QMC energies in 
Fig. 8(a), v ~ 1.3, 1.3, and 2.2 (t fd ~ 0.6, 0.6, and 
0.2) for (n) = 2.0, 2.2, and 2.6, respectively, coincide 
with the onset of rapid increases in magnitude of the sin- 
glet correlator and with similar rapid decreases in the 
low-temperature entropy as discussed earlier. Moreover, 
from the (n) — 2.6 case, we know that all these features 
occur and coincide in the absence of antiferromagnetic 
correlations, so that their fundamental origin is likely as- 
sociated solely with the screening of the local moments 
by the conduction electrons. 
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FIG. 8. (a) The dependence of the QMC and PMHF in- 
ternal energies relative to the AFHF as a function of the in- 
verse-square-root hybridization v 
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Parameters are 



the same as in Fig. 2(a). For all three fillings, the QMC in- 
ternal energy agrees with the AFHF internal energy at large 
v, breaking away at an intermediate v value. This break- 
away, which happens at v ~ 1.3 (t/d ~ 0.6) for (n) = 2.0 
and 2.2, and at v ~ 2.2 (tfd ~ 0.2) for (n) = 2.6, might have 
implications on the volume-collapse phenomenon, (b) The 
FLEX results for the same parameters, with the atomic-limit 
(tfd = 0) offsets removed. The agreement with QMC is quite 
reasonable although the breakaway is more gradual than seen 
for the QMC results in (a) . 



As noted earlier the FLEX approximation does not 
produce the Mott gap and therefore may overestimate 
Kondo-like excitations for smaller tfd, which never dis- 
appear from the spectrum, even at the atomic (strong- 
coupling) limit. This failure of FLEX leads to an under- 
estimation of the internal energy as seen at tfd — in 
Fig. 6, due to the spurious presence of Kondo-like excita- 
tions at this limit. Nevertheless, it is interesting to focus 
only on the tfd dependence and remove these atomic- 
limit offsets by using AE(tf d ) = E(tf d ) - E(0) and plot- 
ting AE — A.Eafhf m the case of FLEX, as is done in 
Fig. 8(b). As may be seen, the agreement with Fig. 8(a) 
is really quite good, although, the breakaway from the 
AFHF is more gradual than in the QMC. It is tempt- 
ing to attribute this more gradual behavior in Fig. 8(b) 
to the persistence of Kondo-like excitations in FLEX to 
the largest v = > whereas presumably such excita- 
tions cease for the QMC results in Fig. 8(a) beyond some 
critical values of v, possibly where the QMC and AFHF 
curves become essentially degenerate. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have extended our previous studies 20 
of the periodic Anderson model (PAM) to nonsymmet- 
ric cases, examining the low-temperature dependence of 
a variety of quantities as a function of the hybridiza- 
tion tfd between the / and d orbitals: local moment, 
singlet correlator, antiferromagnetic structure factor, en- 
tropy, and internal energy. Results from the determinant 
quantum Monte Carlo (QMC), fluctuation-exchange ap- 
proximation (FLEX), and Hartree Fock (HF) methods 
have been compared. 

Beyond exploration of the properties of the nonsym- 
metric PAM itself, these results provide further insight 
into the effect of more exact treatments of correlations on 
mean field approximations to phase transitions driven by 
changes to the ratio of atomic hybridization to Coulomb 
repulsion, such as the volume-collapse transitions in the 
rare-earth metals. Indeed this problem has also moti- 
vated many of our parameter choices for the PAM, which 
were generally selected to maintain approximately one / 
electron per site while varying the number of d electrons 
and consequently the total filling (n). 

One of the fundamental results of this work is the coin- 
cidence of systematic changes in the measured quantities 
across particular values of the hybridization, tfd ~ 0.6, 
0.6, and 0.2 for (n) = 2.0, 2.2, and 2.6, respectively. For 
tfd less than these values, the local moment (~ 1), singlet 
correlator (~ 0), and /-subsystem entropy (~ In 2) are all 
approximately constant, while the internal energy is well 
approximated by AFHF, i.e., a regime of well- formed, un- 
screened local moments. For tfd greater than these val- 
ues, the local moment and entropy drop, the magnitude 
of the singlet correlator grows, and the internal energy 
falls below that of AFHF. The singlet correlator suggests 
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Kondo-like effects in this regime, although, there is a 
complete absence of antiferromagnetic correlations only 
at the largest filling, (n) — 2.6, where this identification 
appears unambiguous. It should be acknowledged that 
changes with tfd for (n) = 2.6 are also somewhat more 
gradual than those closer to half filling, and for this rea- 
son we used slightly different criteria for the crossover t fd 
values: the first significant breakaway from the constant 
values for (n) = 2.6 versus a location of particularly rapid 
change for (n) = 2.0 and 2.2. 

We have presented the correlation energy E — -Epmhf 
most of the way between the strong- and weak-coupling 
limits by varying tfd from to very large values. The 
exact QMC energies interpolate smoothly between the 
AFHF result at the strong- and PMHF result at the 
weak-coupling limit for (n) = 2.6, and we expect simi- 
lar behavior at other fillings in the thermodynamic limit. 
Approximate diagrammatic FLEX agrees very well with 
QMC towards weak coupling and is distinctly improved 
over PMHF by not only reaching the correct weak- 
coupling limit but including the leading order dependence 
on tjl away from this limit. Moreover, the FLEX curve 
is smooth like QMC and generally of the same shape 
throughout the full tfd range, although, it increasingly 
underestimates the internal energy as the strong-coupling 
limit is approached since the FLEX approximation does 
not produce the Mott gap. Even so this is a gradual effect 
and the t fd dependence compares favorably to QMC even 
within the strong- and intermediate-coupling regimes, 
< tfd < 2. In general the FLEX results are in re- 
markably good overall agreement with QMC, suggesting 
this approach as a useful approximation for the study of 
/-electron systems. 

Due to the crudeness of the PAM model, it is diffi- 
cult to make a compelling argument on thermodynamic 
grounds that it exhibits some transition like the volume- 
collapse transition in the rare-earth metals. However, 
mean-field HF calculations for the PAM and for realis- 
tic all-valence-orbital treatments of the rare earths both 
yield transitions of the same fundamental nature, and 
the latter are certainly caricatures of the volume-collapse 
transitions. It thus makes sense to ask how the correlated 
and HF results differ for the PAM and to presume sim- 
ilar differences should be obtained in the realistic case. 

— 1/2 

We find that for decreasing tj d (or, volume) the cor- 
related energies break away from the "localized" AFHF 

— 1/2 

results at values of t fd (or, volume) considerably larger 
than where the "itinerant" PMHF energies cross down 
through the AFHF curve and become stable. Moreover, 
this breakaway coincides with the decrease in the local 
moment and /-subsystem entropy, and the increase in the 
singlet correlator. The coupling of the thermodynamics 
with the singlet correlator behavior is certainly consistent 
with the Kondo volume-collapse scenario, 16 while it is in- 
teresting that the associated decay in the local moment 
is reminiscent of the Mott-transition perspective. 17 
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